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Abstract 

We apply the Golec-Biernat Wiisthoff model in the calculation of vector meson photo- 
and electroproduction. Starting from very simple non-relativistic wave functions we show 
that the model provides a good description of J/^' cross sections in wide and W 
ranges. For the light mesons one obtains the approximately correct W dependence and 
ratio of longitudinal to transverse cross sections, although in this case the normalization, 
affected mainly by the wave function employed, is not in good agreement with data. 

1 Introduction 



We consider elastic vector meson (VM) production in 7*p scattering in the context of the 
Golec-Biernat Wiisthoff (GBW) model |2|. The model is constructed in the proton 
rest frame, as shown in Fig. 1, where the qq pair resulting from a photon fluctuation can 
have a long enough lifetime to allow its interaction with the proton. In this frame the 
qq longitudinal momentum is quite large so that f, the bi-dimensional separation of the 
pair, doesn't change significantly during the scattering. It is then natural to perform the 
calculations in the (r, z) representation. 

The dynamics of the interaction between the qq pair and the proton are absorbed in the 
dipole cross section, aqq-p. This approach to describing the scattering has been followed 
by many authors (see for example @, ^ |5|, ^), and is often called the Color Dipole 
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Figure 1: Scheme for vector meson production in the Color Dipole model. 



model. In the version proposed by K. Golec-Biernat & M. Wiisthoff, the pair-proton 
cross section assumes the form 



(Tqq-p{Q, W,r) = ao\l- exp 



RoiQ^W^] 



ARl{Q^,W^) 



where 



GeV 



1 



4m^ 



A/2 



(1) 



(2) 



One important characteristic of this parameterization is that the cross section does not 
increase indefinitely with r, but "saturates" when r increases beyond Rq. Rq therefore de- 
fines a saturation radius, which depends on Q^, W'^, and the mass of the quark/antiquark. 
The W dependence of this cross section depends on the ratio v/Rq: for values of the ratio 
> 1, the W dependence is approximately flat, while for small values, (yqq-p scales approx- 
imately as VF^^. We therefore expect a steeper T^-dependence for heavy vector mesons 
relative to light vector mesons at a given , W"^ . 

The three free parameters (cio. A, xq) in the expression above were determined via fits 
to DIS (F2) data M. The total 'j*p cross section is calculated as M: 



c7l,t{Q,W) = J d^rj dz|^'L,T(Q,^,r)|V,5_p(g,M^,r) 



(3) 



where ^ l,t{Q^ z, r) stands for the qq component of the wave function of a longitudinally 
(L) or transversely (T) polarized photon. The proton structure function i*2, the data from 
which is used in the extraction of the parameters, is related to (Tl^t via 



47r2 



a 



The numbers resulting from the fit are do = 23.03 mb, A = 0.288 and xq - 
only three light quarks are assumed for the photon wavefunction, and ctq 
= 0.277 and xq = 0.41x10"'^ if one also considers charm. 



(4) 

3.04x10"'' if 
= 29.12 mb, A 
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In these fits the hght quark (u, d, s) mass is used as a parameter and fixed to mu^d,s = 
140 MeV; this average value was necessary to allow the model to give rough agreement with 
the total photoproduction cross section. The charm quark mass used in the calculations 
is discussed below. 

The model has been applied in the calculation of other inclusive processes, and has 
provided good results with no extra parameters. For example, it gives the qualitatively 
correct behavior for the quantity (ii*2 /d In at low Q^. Furthermore the diffractive cross 
section a'^^^ ^ is reasonably reproduced in the model |^ . 

Our aim in this work is to test the model for exclusive processes. We analyze the 
model in photo- and electroproduction of J/^, p and mesons. We first describe the 
calculations, and then compare the results with published HERA data. 



2 Vector Meson Production 

In this section we give the principal formulae for the process 7(*)p — > Vp in the context of 
Color Dipole model. One writes the amplitude for the process os A = <V \ CTqq-p \ 7*^*-* >, 
where |y > is the vector meson wavefunction and I7* > is the qq component of the 7* 
wavefunction. In the (f, z) representation A is given by 

Al,t = J^dzJ d^r ^;^"^(r, z) aqq.p{Q, W, r) ^:^'^(Q, r, z) (5) 

where aqq-p{Q,W,r) is the cross section for the pair-proton interaction. The forward 
scattering cross section is then given by 

da. 



dt lOTT 

Since the model does not incorporate any t dependence we assume the ordinary exponen- 
tial t dependence, da/dt ~ Ce"'^'*' observed in the data, so that the total cross section is 
given by 



The qq component of the photon wave function is quite well known |10|. For photons 
with longitudinal polarization we have 

(27r)3/2 

and for transverse polarization 



i{2z - 1% ^ by \ df{r,z) 



±bx + i{2z - l)by I dr 



h,,-x,}. (9) 
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In the expressions above is the number of colors, eq{mq) is the quark charge (mass), 
Ai^2 are the hehcities of the q and and (p{r, z) is the photon spatial wave function given 
by 

(p{r, z) = Ko{r^Q^z{l-z)+ml ). (10) 

In analogy with the photon wave function one can write for the longitudinally polarized 
vector mesons 

^^ = -2Mv<i)^{z,r)5x,,^x, (11) 
and for the transversely polarized vector mesons 

1 ( i{2z - 1% Tby \ dct)V{r,z) ^^ ^^^^ 



z{l - z) \ ±bx + i{2z - l)hy J dr 

In the case of the vector mesons, the spatial part of the wave function is not known, 
and some hypothesis must be made. In the following sections we discuss the vector meson 
wave functions. 



3 J/^ wave function 



Since the internal motion is slow for heavy quarkonium, we can choose a non-relativistic 
wave function. For transversely polarized photons, Jz = ±1- If we assume that the orbital 
angular momentum of the quark and antiquark in the J/ip is 0, then Sz = ±1. This implies 
Ai = A2, and the second term in Eq. (12) disappears (see also |jll|]). We analyze some 
possible wave functions in {z,kt) space, and then calculate the {z,r) wavefunction using 
a 2-D Fourier transform, 

<f>v{r,z)= J ^ cf>v{z,h)e^^-'^^ . (13) 

3.1 ^-function 

The simplest case one can have is to consider that the q and q have the same longitudinal 
momentum fraction and that the transverse momentum is very small. In terms of a {z, kf) 
space wavefunction, this hypothesis yields 

(l)v{z,kt)=K 5{z-l/2)5Hkt) ■ (14) 

The only free parameter is the normalization, K. It can be determined if we fix the partial 
width for J/ip ^ e~^e~ to the experimentally measured value. 



^=V=^li^^-i^*v(..«l^ (15) 
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resulting in 



K = \\^^ -. 16 

2Mv eqa V 327r ^ ' 



The wave function in the (r, z) space representation, obtained via Eq. (ffsl), is then 



ZMy Cqa V "J^TT 

The full expression for longitudinal and transverse wave functions is obtained by inserting 
the expression above in Eqs. (11) and (12). 

3.2 Gaussian function 

Another possibility is to consider again the same longitudinal momentum fraction for both 
q and q but now we assume that some small transverse momentum is allowed inside the 



J /'if. Following [12 1, we use a Gaussian distribution around zero and the wave function 
reads 

^v{z,kt) = K 6{z- 1/2) eM-a^)- (18) 

We take for a the value obtained from lattice QCD calculations, a = 4.68 and rric =1.43 
GeV. 

Again the normalization is obtained from the decay width. 



ml2Mv eqa V 327r ^ ' 

and the (r, z) representation of the wave function comes from the Fourier transform. 

3.3 Double- Gaussian function 

We propose now that not only the transverse momentum may have some distribution 
but also the longitudinal momentum fraction has a Gaussian distribution around 1/2, 
resulting in 

ct>v{z, h) = K exp[-c2(z - 1/2)2] gxp(-a^). (20) 

Now the wave function has one extra parameter, c. However, one extra constraint is 
available since now the wave function can be normalized. From the decay width 



ml2Mv eqa V 327r ^ erf (c/2) ' ^ ' 
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and from the normalization condition, 

we get (? = 27.2. In this way the (r, z) representation of the wave function is obtained 
via the Fourier transform with no free parameters. 



4 Light meson wave functions 

For the hght mesons, a relativistic approach must be used. We fohow the procedure 
outUned in |]l3|, |14| , where one writes the wave function in terms of the Ughtcone invariant 
variable p"^ (3-momentum of the quark in the non-relativistic limit), 

p' = \{M--^m% ^^ = ^1^- (23) 

In this expression M is the invariant mass of the qq system, and ruq = 140 MeV is the 
quark mass. Again a Gaussian form is applied, now in the p variable, 

1 



(l){p) = K (2^i?2)3/2 



(24) 



The value of R was fixed by requiring that the exponential gives a value of 1/e for M = 
My. This yielded 

Rl = 15.5 GeV-2, 
Rl = 8.3 GeV-2. 

Using Eqs. (23) above one obtains the wave function in the {z,kt) representation, and 
through Eq. (13) 

<A(r, z)=K 4z(l - z) V2^exp[-^-^^] exp[ ^ M exp[ 1-] . (25) 

The normalization factor K is determined from the normalization condition, so that once 
again there are no free parameters. 



5 Results 

We have fixed the u, s quark masses to 140 MeV, as in [|l|, and use either the three or 
four quark parameters given in the introduction (for the J/V' case, only the four quark 
parameters are used). For the slope B (see Eq. (|^)) we use a parameterization to experi- 
mental data, given by ||l5| 

B = 0-»(^g.^^ + l) . (20) 
This parameterization is used for all the vector mesons. 
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5.1 J/^ production 



In Fig. 2, the W dependence of the 7*p cross section for J/'^ elastic production is shown 
for various values of . It can be seen that even for the simplest wave function utilized 
(5-function, dotted line) the W dependence predicted by the model is in agreement with 
the data, although for photoproduction the normalization is not very good for this wave 
function. For the more complicated forms we obtain good agreement to HERA data both 
in shape and in normalization. Note that the value of Rq in Eq. (§) is never very small in 
J /tp production because of the presence of the charm quark mass, implying that we are 
far from the low or high W saturation regions. The W dependence is then given to 
good approximation by cr^*p^j/ipp oc W^^. 




Figure 2: The W-dependence of elastic J/\I/ production as a function of . Open circles 
are ZEUS data P^ , [7^ and solid circles HI data / |7^ The dotted lines correspond to 

calculations using the 6-function wave function, the dashed lines to calculations using a single 
Gaussian, and the solid lines to the double Gaussian wavefunction as discussed in the text. The 
values, in units of GeV^, are indicated in the figure. The data and calculations have been 
scaled by the values given for clarity of presentation. 
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Fig. 3 gives the ratio of the longitudinal to transverse cross sections, (Ji/aT, for fixed 
W as a, function of (upper plot), and for fixed as a function of W (lower plot). Due 
to the delta function in z for the 6 and Gaussian wave functions, the integration over cfrdz 
is the same for transverse and longitudinal amplitudes, up to a constant factor. In this 
case the model dependence for R disappears and one gets just R = Q2M2/(l6m^). The 
dotted and dashed curves are therefore superimposed in Fig. 3. For the double Gaussian 
case this does not occur; Al is peaked near z = 1/2 and the growth is damped relative 
to At- It is clear that the dependence of i? = (Jl/cft is quite sensitive to the vector 
meson wavefunction, while there is little sensitivity in the W dependence. The data favor 
the double Gaussian wavefunction. 




Figure 3: Ratio R = (Jl/cjt for elastic J/^ electroproduction. The dotted and dashed lines 
(superimposed) correspond to the 6 and single Gaussian wave functions. The solid lines are 
calculated with Double Gaussian wave function. The upper plot gives the dependence at 
fixed W = 90 GeV. The lower plot gives results versus W for = 6.2 GeV^ (upper curves), 
and = 0.^5 GeV^ (lower curves). The open circles are ZEUS data / [7^/ and the solid circles 
are HI data 
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5.2 p production 

The comparisons to data for elastic p production are given in Figs. 4-6. The hghter mass 
of the p allows the cross sections to enter the saturation region (r ^ Rq), and the W 
dependence is a strong function of Q"^, as seen in Fig. 4. The two sets of curves on the 
graph give the results for three quark (dashed) and four quark (solid) parameter sets. The 
normalization of the cross sections is not particularly good, and the dependence is also 
not well reproduced, although the 3 parameter set does a better job of reproducing the 
dependence. Note that the normalization is affected both by the choice of p wavefunction 
and by the assumed B values. However, the observed change of the W dependence is 
rather well reproduced. This is quantified in Fig. 5, where the power 5 from a fit of the 
W dependence with a form a^*p^pp oc is compared with data. 

The dependence oi R = gl/cft on at fixed W is shown in Fig. 6, as is the W- 
dependence at fixed Q^. The model, with the simple wavefunction chosen, gives a good 
representation of the data. In this case, the model predicts a small variation of R with 
W. 

5.3 (j) production 

For the ^ meson, we compare the predicted W dependence for 7*p — > as a function of 
with data in Fig. 7 and R versus at fixed W, as well as R versus W at fixed in 
Fig. 8. The normalization of the cross section is too high given the chosen wavefunction 
and parameterization for but the R dependence on is in reasonable agreement with 
the available data. 

6 Summary and Conclusions 

We have used the Golec-Biernat Wiisthoff model to calculate the cross section for exclusive 
vector meson production in ep scattering. The model gives the correct W dependence for 
both photo- and electroproduction of J/tp mesons. This feature can be attributed to the 
model because it is observed even when a very simple wave function is employed. This was 
not unexpected, since the model reproduces the pQCD expectations for the Q"^ and W 
dependence of a^*p^^p when the dipolc size is small compared to the "saturation radius" . 
Perhaps more surprising is the fact that a simple J/ip wave function (double Gaussian) 
also gives good results for the normalization of the wavefunction as well as for the 
dependence oi R = ul/c^t- 

For p and (f) production, the W dependence changes dramatically as we go from photo- 
production to DIS. The change in the W dependence versus calculated using a simple 
wavefunction can reproduce the data reasonably well. The behavior of R versus and 
W is also in reasonable agreement with the data. However, the normalization of the cross 
section is not in agreement with the data, and the discrepancy with data varies with Q^. 
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This problem can be attributed to the lack of knowledge of the wave function and perhaps 
also in part to the change in t-slope with (this slope is needed to calculate the total 
cross section). 

The proton rest frame description of ep scattering gives a simple and intuitive picture 
of low and small-x scattering. Golec-Bicrnat and Wiisthoff have shown, using a three 
parameter parametrization for the qq—p scattering cross section, that the total as well as 
the difFractive cross sections can be well described. We find that the the model can also 
describe many of the features of exclusive vector meson production. 
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Figure 4: The W- dependence of elastic p production as a function of compared with recent 
ZEUS [7^/ and HI p^ , [^/ data. The values are given in the plot. The solid lines are 
for the four quark parameters, while the dashed lines are for the three quark parameter sets. 
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Figure 5: The W -dependence of elastic p production as a function of Q"^, compared with recent 
ZEUS / |7^ (open circles) and HI (solid circles) data. The solid line is the result using the 
four flavor parameter set, and the dashed curve is the result using the three flavor parameter 
set. 
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Figure 6: R dependence on and W in elastic p production. The upper plot gives the 
dependence at fixed W = 75 GeV, while the lower plot gives the W dependence at the values 
indicated in the plot. The open circles are ZEUS data / |7^ and the solid circles HI data f^. 
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Figure 7: a^*p^^p as a function of W for different Q'^ , compared with ZEUS data ^^] . The 
solid lines are for the four flavor parameter set, while the dashed lines are for the three flavor 
parameter set. 
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Figure 8: Ratio of longitudinal over transverse cross sections in production. The upper plot 
shows the dependence at fixed W , while the lower plot shows the W dependence at fixed . 
The solid lines are for the four flavor parameter set, while the dashed lines are for the three 
flavor parameter set. The data are from ZEUS jl^/ and HI ^djj . 
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